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Propagators of diagonal and off-diagonal gluons are studied numerically in the Maximally Abelian gauge of 
the SU(2) lattice gauge theory. We have found the strong enhancement of the diagonal gluon in the infrared 
region. The enhancement factor is about 50 at the smallest available momentum, 325 MeV. We discuss also 
various analytical fits to the propagators. 



Propagators of fundamental fields are basic 
quantities in quantum field theories. In non- 
Abelian gauge theories the infrared behaviour of 
the gluon propagators is believed to carry an in- 
formation about the color confinement. A popu- 
lar approach to the confinement is based on the 
dual superconductor mechanism T* which is most 
thoroughly studied in the Maximally Abelian 
(MA) gauge |2] (see, e.g., [3] for a review). A 
study of the MA gauge gluon propagator in the 
coordinate space j4j has shown that the propaga- 
tor of the off-diagonal gluons at large distances 
is exponentially suppressed with respect to the 
diagonal propagator by the effective mass about 
1.2 GeV. This property of the gluon propagators - 
supporting the Abelian dominance in gluodynam- 
ics |5I6| - was confirmed and further explored in 
the momentum space in Refs. |7I8) . The mass gap 
generation for the off-diagonal gluons was also in- 
vestigated analytically ||9j. Below we provide nu- 
merical results supporting the Abelian dominance 
in terms of the gluon propagators. 

We study the SU(2) pure gauge model with the 
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Wilson action. We use the standard parameteri- 
zation of SU{2) link matrices Uu = cos ipe^^ and 
Ui2 — sin ipe^^. The gauge fields are given by 

A'^^{x)a^ = -z{U,,{x)-Ulix)). 

where cr" are the Pauli matrices. In terms of the 
link angles one gets^: 

^lix) = 2 sin (/J^ (a;) sin X,, (a;) , 

^li^) = 2sin(^p(a;)cosxp(a;) , (1) 

Al{x) = 2cos(p^(a;) sin6'^(a;) . 

We call A'^{x) the diagonal gluon field, and 
A^^{x), i = 1,2, the off-diagonal gluon field. 
The MA gauge is fixed by maximizing 

^mag[^]=E cos 2^^ (a;). (2) 

This condition fixes the SU(2) gauge group up 
to a U(l) subgroup. In order to fix the remain- 
ing U(l) gauge symmetry we implement a gen- 
eralization of the Landau gauge maximizing the 
functional 

^La„*d ^] = E,.^ '^"^ (^) > (3) 

which is consistent with the definition of the 
field. In the continuum limit the definition JSJ 
corresponds to the standard U(l) Landau gauge. 

^Note that in Ref. [7| the definition of the field A differs 
from Eq. Q by the factor of 2. 
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We calculate the diagonal and ofF-diagonal 
propagators, 

D-^^^Hp) = {Al{k)AU-k)) , (4) 
D°f'^Hp) = {A+{k)A-i~k)), (5) 

respectively. Here — {Aj^ ± «74^)/V2 and the 
Fourier transformed field, A^^{k), is defined as 

A!^{k) = e-^'^^---^'^- . (6) 

The lattice momentum variable is — 
(2/a) sin(aA;^/2), where fc^ = 2'Kn^/{aL^) with 
= 0, — 1, /X = 1, . . . , 4. In terms of 
p the lattice propagator of a free massive scalar 
particle in momentum space has a familiar form, 
D{p) (X l/lp^+m?). The local gauge condi- 
tion corresponding to the maximization of Eq. Q 
takes a simple form, p^A^^ — 0. 

The most general structure of both diagonal 
and off-diagonal propagators is 

M = (5.. - ^) Ab^) + ^ A(P^) , (7) 

where Dt i are the scalar functions. We have three 
independent formfactors, Df"*® and be- 
cause in the Landau gauge Dj'^*"" = 0. 

Below we present the results on 32'* lattice at 
(5 — 2.40 which corresponds to the lattice spacing 
a — (1.66 GeV)^^. We generated 30 configura- 
tions which are fixed to the MA gauge by the 
Simulated Annealing algorithm |10) with 10 ran- 
domly generated gauge copies. Then the Landau 
gauge is fixed by the local over-relaxation algo- 
rithm with 20 randomly generated gauge copies. 

We show the gluon dressing functions, 
p^ Dt^i{p^), vs. momentum, p, in Figure ^ The 
diagonal (transverse) dressing function has a rel- 
atively narrow maximum at non-zero momentum 
Pq^'^^ « 0.7 GeV. Its behavior at small mo- 
menta is qualitatively very similar to the behav- 
ior of the gluon propagator in the Landau gauge 
(see, e.g., JI]). The longitudinal part of the off- 
diagonal dressing function has a wide maximum 
at pg^'^"'^ « 2 GeV. The transverse off-diagonal 
dressing function is a monotonically rising func- 
tion for all available momenta. 

The formfactors Df'^'''^{p^) and of^'^^ip^) 
coincide at small momenta. This implies that in 



the IR region the off-diagonal propagator is 
Dff^Hp) « S.^-Df^^^ip') , / < 1 GeV, 
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* Transverse diagonal 

■ Transverse off-diagonal 

• Longitudinal off-diagonal 



Figure 1. The gluon dressing functions p^ D{p^ 
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Figure 2. The diagonal/off-diagonal ratio ©. 

The diagonal formfactor is dominating over the 
off-diagonal ones. In Figure HI we plot the ratio 

i?(p2)^7^f-g(p2)/^offdiag(^2)^ (8-) 

which reaches the value 50 at the smallest avail- 
able momentum, p = 325 MeV, suggesting that 
the diagonal gluons are responsible for infrared 
physics. At higher momenta the suppression of 
the off-diagonal propagator becomes weaker. 

We fit the formfactors at low momenta by the 
functions 



(fit 1) 



Dip') 
D{p') 



p2{l + ^)^^2{l + c.) •> (fit 2) , 



(9) 
(10) 

(Yukawa fit) , (11) 



Dip")^ , (Yukawa 2 fit), (12) 

Dip")^ (Gribovfit), (13) 
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where Z, a , m and n are fitting parameters. 



diagonal, transverse 




Figure 3. Fits P <10I12I13() of the transverse di- 
agonal formfactor at low momenta. 



fit 


TO, GeV 


a or K 


Z 


fit 1 


0.73(2) 


0.92(3) 


16.9(4) 


fit 2 


0.58(2) 


0.49(5) 


8.5(2) 


Gribov 


0.33(1) 




4.58(5) 


Yukawa 2 


0.50(2) 


0.19(3) 


8.3(3) 



Table 1 

The best parameters of the low momentum 
fits of the diagonal propagator. 

Here we concentrate on the fits of the diago- 
nal propagator, shown in Figure |3| The best fit 
parameters are presented in Table ^ The three 
parameter fits H9I10I12|I are working well in a low 
momentum region (p,„ax 0.8,0.4,0.9 GeV, re- 
spectively). The corresponding curves are almost 
indistinguishable from each other. The mass pa- 
rameter TO do not coincide in these fits, and the 
difference between its values is about 30%. We 
have also applied the two-parameter fits given 
by Yukawa Hll|) and Gribov (|13|l functions. The 
Gribov fit is working well in the region Pmax ^ 
0.9 GeV. The Yukawa fit of the diagonal propa- 
gator does not work at all (we get jd.o.f. ~ 6 
for fits in Pmax < 1 GeV region). 

Summarizing, our results obtained in the MA 
gauge of SU{2) gluodynamics clearly show that 
in the infrared region the propagator of the di- 
agonal gluon is much larger (about 50 times at 
lowest available momentum) than the propaga- 
tor of the off-diagonal gluon. Thus the colored 
objects at large distances interact mainly due to 



an exchange by the diagonal gluons in agreement 
with the Abelian dominance |5lf)| . The diagonal 
propagator has qualitatively the same momentum 
dependence as the gluon propagator in Landau 
gauge. All infrared fits suggest that the diagonal 
propagator contain massive parameters, although 
it is not of the Yukawa form. Further results on 
the propagators in the MA gauge can be found in 
Ref. [g. 
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